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We have performed a search for the flavor-changing neutral-current decays B — » tt£ + £~, where 
£ + £~ is either e + e~ or fi + fx~ , using a sample of 230 x 10 6 T(4S) — > B-B decays collected with 
the BABAR detector. We observe no evidence of a signal and measure the upper limit on the 
isospin-averaged branching fraction to be B(B — > tt£ + £~) < 9.1 x 10 -8 at 90% confidence level. We 
also search for the lepton-flavor-violating decays B — > Tre*/!^ and measure an upper limit on the 
isospin-averaged branching fraction of B(B — -> ire fj,^) < 9.2 x 10~ 8 at 90% confidence level. 

PACS numbers: 13.20 He 



In the Standard Model (SM), the decays B -> Tr£ + i~, 
where £ + £~ is either e + e~ or proceed through 6 — ► 

d£ + £~ flavor-changing neutral-current processes (FCNC) 
that do not occur at tree level. Three amplitudes con- 
tribute at leading order: a photon penguin, a Z penguin, 
and a W + W~ box diagram. With highly suppressed SM 
rates, predicted to be (3.3 ±1.0) x 10 -8 these decays 
provide a promising means to search for effects of new 
flavor-changing interactions. Such effects are predicted 
in a wide variety of models, usually in the context of 
b -> s£+£- decays @, S The b -> d£+£~ decay in- 
volves quark-flavor transitions different from b — > s£ + £~ 
and thus its measurement constitutes an independent 
test for new flavor-changing interactions. An experimen- 
tally similar but otherwise unrelated process, the lepton- 
flavor-violating (LFV) decay B — > Tre ± /i T , is forbidden 
in the SM, but can occur in some models beyond the 
SM, such as theories involving leptoquarks Earlier 
searches by other experiments [f| have reached branch- 
ing fraction upper limits at the 10 -3 level for the FCNC 
decay and the 10~ 6 level for the LFV decays. 

In this Letter we report the findings of a search for 
the decays B -> nf+fr and B -> ire ± ^ in 208.9 ftT 1 
of data recorded at the T(AS) resonance, corresponding 
to (230.1 ± 2.5) x 10 6 BB decays. The data were col- 
lected with the BABAR detector [6j at the PEP-II storage 
ring located at the Stanford Linear Accelerator Center. 
The event selection criteria are optimized using simulated 
data and data samples independent of those selected as 
signal. The signal model used for efficiency evaluation 
of the £ + £~ modes uses form-factors from [7| and ampli- 
tudes from Calculations of the same type have pre- 
viously been shown by Q to describe the kinematic dis- 
tributions of B — > K£ + £~ well. The efficiency of the efx 
event selection is estimated using a 3-body phase space 
model with QED photon radiative corrections. 

We reconstruct signal events by combining two oppo- 
sitely charged leptons (e + e~, f-i + fi~ or e^fjF) with a pion 
(tt^ or 7T°). Electron (muon) candidates are required 
to have a momentum larger than 0.3 (0.7) GeV/c. We 
suppress backgrounds due to photon conversions in the 
B — > 7T e + e~ channels by removing e + e~ pairs with in- 
variant mass less than 30 MeV/c 2 . Bremsstrahlung pho- 
tons from electrons are recovered if the photon has an 
energy of E > 30 MeV and a direction within a small 
angular region around the initial electron momentum vec- 
tor. The identification of electrons (muons) is about 92% 



(68%) efficient on average with a hadron misidentification 
rate of less than 1% (4%). Charged pion identification 
is more than 85% efficient and has a kaon misidentifica- 
tion rate of less than 5%. Neutral pions are identified 
as pairs of photons, each having an energy of at least 
50 MeV. The invariant mass of the pair is required to 
satisfy 115 < ra 77 < 150MeV/c 2 . 

Correctly reconstructed B decays produce narrow 
peaks in the distributions of two kinematic variables: the 
beam-energy substituted mass, tties = ^/i?^ 2 — \p* B \ 2 , 
and AE = E* B — E£. Here, E£ is the beam energy 
and E B (p* B ) is the energy (momentum) of the recon- 
structed B meson, evaluated in the center-of-mass (cm.) 
frame. For signal events the ibes distribution is cen- 
tered at the B-meson mass and the AE distribution is 
centered at zero. The mean and width of these distri- 
butions are determined from smearing and shifting the 
values from simulated signal events according to studies 
of B + — > J/ipK + and B° — » J/tpir events in data control 
samples and simulations. We find the width of Toes to 
be 2.5 (1.8) MeV/c 2 for the tt ± (vr°) modes and widths 
of AE to be 23, 50, 20 and 39 MeV for the % ± e + e~, 
7r°e + e~, 7r ± ^ + /i~ and 7r°/i + ^~ final states, respectively. 
For events reconstructed as e^/x^, we assume the same 
mean and width as for the corresponding e + e~ modes. 

The primary sources of background are random com- 
binations of particles from e + e~ —> qq (q = u, d, s, c) and 
from Y(4S) — > BB decays. These combinatorial back- 
grounds typically arise from pairs of semileptonic decays 
of B and £>(*) mesons. Additionally, there is background 
from events that are peaking in /ties and AE as they 
have the same topology as signal events. These events 
include B -► J/rp-K {J/ip -> B ± -> J/i(;K ± or 

B ± -> (with K ± misidentified as tt ± ), and 

B — > irhh (with two hadrons h = K , ^ misidentified 
as muons). 

Contributions from e + e _ — > qq processes are reduced 
by exploiting the difference between the spherical track 
distribution in BB events and the jetlike structure of 
e + e~ — > qq events. We consider events for which the 
ratio of second to zeroth Fox- Wolfram moments R2 is 
less than 0.5. Further suppression by a factor of ~ 45 
is obtained by constructing a Fisher discriminant from 
the following four quantities [ij defined in the center-of- 
mass frame: R%, \ cos 9thr \ where 8thr is the angle between 
the thrust axis of the signal particles and that of the 
remaining particles in the event, | cos 8 B \ where 8 B is the 
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angle of the B candidate's momentum vector with respect 
to the beam axis, and the ratio of second- to zeroth-order 
Legendre moments 10 . 

Combinatorial background from BB events is reduced 
by a factor of ~ 3 by using a likelihood ratio composed 
of : the missing energy of the event (computed from all 
charged tracks and neutral energy clusters), the vertex fit 
probability of all tracks from the B candidate, the vertex 
fit probability of the two leptons, and cos 8b- Missing 
energy provides the strongest suppression of these events, 
which typically contain energetic neutrinos from at least 
two semileptonic B or meson decays. 

We veto events that have a dilepton invariant mass 
consistent with the J/ip resonance (2.90 < m e + e - < 
3.20GeV/c 2 and 3.00 < m p+A1 - < 3.20GeV/c 2 ) or with 
the ip(2S) resonance (3.60 < m t + e - < 3.75GeV/c 2 ). For 
electron modes, the vetoes are applied to mi+i- com- 
puted both with and without bremsstrahlung recovery. 
When a lepton radiates or is mismeasured, mi+i- may 
shift to values below the charmonium mass, with AE 
shifting downward accordingly. Therefore, we veto events 
that lie in linearly dependent AE-mu bands, whose 
widths are determined from simulation, similar to the 
technique applied in Q. For e ± /i T modes, we use the 
same vetoes as for the e + e~ modes. In modes with 
muons, in order to veto events with tracks that are con- 
sistent with hadronic decays D — > Kir or D — * tttt, we 
require mu and m^e to lie outside the range 1.84 — 1.89 
GeV/c 2 when the £ is assigned the mass of a n or K. For 
the 7T° modes, the range for m^i is increased to 1.79 — 1.94 
GeV/c 2 . 

The events removed by the charmonium vetoes are 
kinematically similar to signal events and serve as large 
control samples for studying signal shapes, selection ef- 
ficiencies, and systematic errors. The branching frac- 
tion of B — > J/tpn is also extracted from the control 
sample and found to be in agreement with the current 
world average [llj. We also select a control sample of 
B + — > J/ipK + events to measure the efficiencies and 
systematic uncertainties of lepton identification and the 
Fisher and likelihood selection. 

We extract the signal yield by counting events within 
a signal region defined as ±2cr around the mean values of 
the Toes and AE distributions expected for signal events, 
and comparing observed event counts with estimations of 
the remaining background in the same region, summa- 
rized in Table HI 

To determine the peaking background from hadronic 
B — > 7T7r7r or B — > Kirn events, we select a control sample 
where one track is required to pass hadron identification 
in place of muon identification. This selects hadronic B- 
decays where the remaining track that passes the muon 
selection is a mis-identified hadron. Each event is fur- 
ther weighted by the probability that one more hadron 
is mis-identified as a muon. The expected contribution 
to the B — > tt£ + £" signal region is extracted from a one- 



dimensional fit to the distribution of toes for events that 
pass the AE signal selection. 

Backgrounds from B decays to final states with real 
leptons are estimated from high-statistics samples of sim- 
ulated events. B — » K£ + £~ events are the largest peak- 
ing background component for the 7r + e + e~ mode, but 
are shifted towards lower AE than signal and fall out- 
side the signal region. B — > p£ + £~ events contribute 
even less, since the reconstructed B mesons are missing 
a pion. Background from charmonium resonances are 
found to be negligible. 

The expected number of combinatorial background 
events is extracted from a two-dimensional, unbinned 
maximum-likelihood fit to toes and AE in a sideband 
defined by 5.2 < m ES < 5.2724 GeV/c 2 and \AE\ < 
0.25 GeV, i.e., below the toes value expected for signal 
B events. The signal-region yield is obtained from ex- 
trapolation of this fit into the signal region. This pro- 
cedure has been validated by studies of simulated back- 
ground events and data events in the efj, channel where 
no signal-like events are expected. The background prob- 
ability distribution function (PDF) is modeled as the 
product of an ARGUS function 12] for toes an d an ex- 
ponential function for AE. The slopes and normaliza- 
tion are floating in the fit. Average biases in the back- 
ground central value and its uncertainty were corrected 
for, based on a study of a large ensemble of simulated ex- 
periments generated from the background PDF obtained 
from data. The corrections amount to 35% in the low- 
statistics B° — > 7r°/i + /i~ channel and < 10% in all others. 

Systematic uncertainties due to the background esti- 
mates are summarized in Table |U The uncertainty in 
the combinatorial-background estimate is determined by 
varying the fit parameters by ±ltr of the best fit. We also 
consider the effect of using alternative PDF parameter- 
izations on the background estimates, and use the com- 
puted differences to bound the systematic uncertainty. 
Alternatives considered include a PDF that is correlated 
in toes and AE via a linear AE dependence in the toes 
slope parameter, and PDFs for which the AE shape is a 
linear or quadratic polynomial. For peaking background 
with real leptons the uncertainty is dominated by lim- 
ited knowledge llj of the branching fractions for these 
processes, and for hadronic B peaking background the 
uncertainty is dominated by the control sample statistics 
from which it is derived. 

Systematic uncertainties due to the signal efficiency in- 
clude: charged-particle tracking efficiency (0.8% per lep- 
ton, 1.4% per charged hadron) and identification (0.7% 
per electron pair, 1.9% per muon pair, 0.5% per pion), 
neutral pion efficiency (3%), the Fisher and likelihood 
selection (1.4% for all modes involving electrons, 1.7% 
for B+ -> n+fi+p- and 1.9% for B° -> ttV^Z-O, and 
signal simulation statistics (0.1%). A systematic un- 
certainty in signal-region selection efficiency arises from 
the uncertainty in the mean and width of the toes an d 



TABLE I: Number of background events with associated systematic uncertainties expected in the signal region. 



7T + fl + fj, 



ttVV 



7r + e/i 



7r ,J e/i 



mES-AE 1 fit 
mES-A_E correlations 
AiS shape 
Peaking (£+£") 
Peaking (hadronic) 

Total 



0.84 ±0.24 
±0.02 
±0.02 

0.057±0.0f6 
< O.OOf 

0.90 ±0.24 



0.43 ±0.23 
±0.03 
±0.0f 

0.009±0.003 
< O.OOf 

0.44 ±0.23 



0.90 ±0.25 
±0.06 
±0.f2 
0.032±0.008 
0.027±0.033 



0.23 ±0.20 
±0.03 
±0.02 
0.005±0.00f 
0.035±0.022 



f.55±0.34 
±0.f7 
±0.3f 
±0.00f 
±0.02 



0.0 
0.0 



f.22±0.43 
±0.05 
±0.24 
±0.00f 
±0.02 



0.0 
0.0 



0.96 ±0.29 0.27 ±0.20 



f.55±0.49 



f.22±0.50 



AE distributions determined from charmonium control 
samples. This contributes a total uncertainty of 0.7% 
for charged modes for which a high-statistics sample of 
B + — > J/ipK + events is used, and a total of 7% un- 
certainty for neutral modes for which a small statistics 
sample of B° — ► J/tpir events is used. For the elec- 
tron modes, we vary the amount of the bremsstrahlung 
tail in the AE distribution, introducing a systematic un- 
certainty of 1-1.4%. The number of BB events in the 
data sample is known to a precision of 1.1%. Addi- 
tional systematic uncertainties for the efficiency result 
from the choice of the form factor model and the relative 
magnitudes of the b — > d£ + £~ amplitudes, which affect 
the distribution of four-momentum transfer q 2 — rnj +£ ^ 
of the signal. We evaluate these systematics from the 
spread in efficiencies when using alternative form-factor 
models 13[ , and when varying the Wilson coefficients in 
the amplitudes by a factor of ±2. The former uncer- 
tainty varies from 1.1% for B + — » 7r + e + e _ to 7.3% for 



B° 
for B° 



7rW 



7T°/i 1 fl 



"; the latter uncertainty varies from 0.3% 



to 1.2% for 



For the 



e/x modes we use the spread in efficiency when apply- 
ing two alternative theoretical models for these decays, 
which amounts to 17% (19%) for the tt ±( ^ mode. The to- 
tal systematic uncertainty of the signal efficiencies are 4% 
(9%), 6% (11%) and 17% (21%) for 7r ± ^°^ee, ir 
and Tr^^e/i modes, respectively. 

Figure [T] shows the distribution of events from data in 
the m-Es~AE plane. The rectangles in the plots indicate 
the signal regions. Three B — > ir£ + £~ candidates and 
one B — > ne^fiT candidate are observed in the signal re- 
gions, which is consistent with the expected background. 
In Table [TT| we calculate the branching fraction upper 
limits at 90% confidence level (C.L.) using a frequentist 
method that takes systematic uncertainties and their cor- 
relations into account. We follow the algorithm of [Tl |. 
but differ from it in that we assume Gaussian distribu- 
tion truncated at zero for the systematic uncertainties 
in signal sensitivity and background expectation. We 
combine modes and determine the e-/i-averaged branch- 
ing fractions to be B{B+ -> tt+£+£~) < 1.2 x 10~ 7 
and B{B° -> ir°£+£~) < 1.2 x lO" 7 at 90% C.L., 
where charged conjugate modes are implied. Defining the 
isospin averaged branching fraction B(B — > Tt£ + £~) = 



B(B+ -> tt+£+£-) = 2x ^B{B° ir a £ + £~), where 

the different B-meson lifetimes tb [H| are taken into ac- 
count, we find the combined upper limit 

B(B -► tt£+£-) < 9.1 x 10~ 8 at 90% C.L. 

This is about a factor three above the nominal SM pre- 
diction [l| . We similarly compute the combined limit for 
the en modes of 

B(B -> Tre ± ^) < 9.2 x 10~ 8 at 90% C.L. 
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FIG. 1: m,ES-AE distributions for events selected in each 
mode. The rectangles indicate the signal regions. 



In conclusion, we have presented the result of a search 
for B -> ir£+£- using a sample of (230.1 ± 2.5) x 10 6 
BB pairs produced at the T(4S) resonance. No excess 
of events is observed in the signal regions, and at 90% 
confidence limit we measure the upper limit of B(B — ► 
ir£ + £~) < 9.1 x 10~ 8 , which is within a factor three of 
SM expectations. We also measure the upper limit of the 
lepton- flavor-violating branching fractions to be B(B — > 
Tre^T) < 9.2 x 10~ 8 . 
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TABLE II: The observed yields, number of expected back- 
ground events, signal efficiency, and branching fraction (B) 
upper limit (U.L.) at 90% C.L. in units of 10 -7 . The upper 
limits for combined modes are also given. 







Observed 


Expected 


Signal 


B U.L. 




Mode 


events 


background 


efficiency 


90% C.L. 


B+ 


— > ir + e + e~ 


1 


0.90 ±0.24 


(7.1 ± 0.3)% 


1.8 


B° 


-> Tv°e+e- 





0.44 ±0.23 


(5.7 ±0.5)% 


1.4 


B+ 


— > n + fj, + n~ 


1 


0.96 ±0.29 


(4.7 ±0.3)% 


2.8 


B° 


-» 7V°ll+fl- 


1 


0.27 ±0.20 


(3.1 ±0.3)% 


5.1 


B+ 


— » ir + e ± fj^ : 


1 


1.55 ±0.49 


(6.3 ±1.1)% 


1.7 


B° 







1.22 ±0.50 


(3.7 ±0.8)% 


1.4 


B+ 


-» TT+£+r 








1.2 


B° 










1.2 


B - 










0.91 


B -> 7re ± /i =F 








0.92 
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